DNA-based tunable THz oscillator 
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. Abstract 

S 

' The intrinsic double helix conformation of the DNA strands is known to be the 

Q I key ingredient of control of the electric current through the DNA by the perpendic- 
ular (gate) electric field. We show theoretically that Bloch oscillations in the DNA 

■ are also strongly affected by such lateral field; the oscillation frequency splits into 

cn ' 

■ a manifold of several generally non-commensurate frequencies leading to a compli- 
Q\ • cated pattern of the charge motion. The frequency of the oscillations falls in the THz 

|: ._,,™,,.-...,.„.,..„,...„.,™._™, 
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1 Introduction 



The helical symmetry of the DNA is almost always neglected when modeling the charge trans- 
port through the DNA-based devices: a DNA is usually considered as a flat ladder-like sequence 
(periodic or stochastic) of base pairs [1,2,3,4,5]. Very recently the intrinsic helix conformation 
of the DNA strands was put forward as the key ingredient allowing for control of the electric 
current through the molecule by the electric field perpendicular to the DNA axis [6]. In the 
presence of such a field the helical conformation leads to an additional periodic modulation of 
the base pair energies, which affects strongly the charge transport through the DNA. On this 
basis, prototypes of the single-DNA-based field effect transistor and the Esaki diode analogue 
have been proposed [6]. 

In this contribution, we aim to further exploit the symmetry of the DNA molecules and demon- 
strate, by means of numerical simulations, that (i) - a single dry periodic DNA molecule (such 
as the ploy(G)-poly(C)) subjected to a coUinear uniform electric field can exhibit Bloch oscil- 
lations [8] and (ii) - because of the helical conformation, the oscillations become more complex 
in the presence of the perpendicular (gating) electric field. The frequency of these oscillations 
falls in the THz domain, the region which is in the focus of an intense research nowadays (see 
Ref. [7] for a recent overview). The above-mentioned property provides therefore for a possibility 
to design a DNA-based tunable THz oscillator. 



2 Model 



Bloch [9] and Zener [10] argued on the theoretical grounds that an electron, moving in an ideal 
periodic potential and subjected to a uniform electric field F, is confined within a finite region 
because of the Bragg reflection. Due to the conflnement, it undergoes a periodic motion which 
is characterized by the angular frequency cub — eFa/h and a spatial extension Lb = W/{eF), 
where — e is the electron charge, F is the applied electric field, a is the lattice constant, and W 
stands for the band width (see also Ref. [11]). 



For simplicity, we consider here a single-stranded uniform poly(G) helix for which we assume the 
conformation parameters of the double-stranded DNA (dsDNA) in its B form, in particular, the 
full- twist period of 10 base molecules. For a periodic dsDNA, such as the poly(G)-poly(C), the 
picture of Bloch oscillations is more complicated as compared to the considered case, however, 
the physics is very similar. We use the minimum tight-binding model, extending it to include 
a uniform electric field E that is tilted by the angle 9 with respect to the axis of the hehx. 
The field has therefore both the coUinear component, F = Ecos9, and the lateral (gating) one, 
£^sin^. The Hamiltonian of our model in site representation reads: 

N N-1 

H^Yl \^n\n){n\ + J ^ |n + + h.c. . (1) 

n=l n=l 

Here, \n) is the state vector of the n-th base molecule and the corresponding energy is given 

by 

TT^H (2) 

where e!^^ is the site energy of the n-th base molecule at zero field. We assume all molecules 
to be the same, so we set the unperturbed energies e!^^ to zero from now on. The term U\\n 
describes the linear potential along the helix axes; the potential drop across a base molecule 
in the stacking direction is f/y = eaE cos 6, where a is the nearest-neighbor distance along the 
helix axis, U± — eEr sin 6 is the potential drop across the helix in the perpendicular direction, 
r ~ 1 nm being the helix radius. The phase which determines the azimuth of the strand 
with respect to the field, is set to zero. The term J in Eq. (1) describes the transfer interaction 
between the nearest-neighbor bases; it is chosen to be positive which implies that the considered 
charge is a hole [12]. 

The parallel component of the electric field, E cos 6, yields the potential ramp along the stacking 
direction, which sets the frequency of the Bloch oscillations, as in the traditional case. However, 
Eq. (2) shows that the helix in the lateral field acquires the additional periodic modulation of the 
potential. This modulation leads to the modification of the electronic structure of the system: 
the bare energy band splits into several different minibands, which is crucial for the charge 
transport properties [6] . Each such miniband has its own Bloch frequency, resulting in a more 
complex overall picture of Bloch oscillations as we show below. The amplitude of the periodic 



modulation is controlled by the magnitude of the perpendicular component of the electric field, 
£^sin^, providing for a mechanism to alter the fundamental properties of the system. 

We further solve the time-dependent Schrodinger equation (the Planck constant h — 1) 

ilpn = SnV'n + J ("^n+l + "^n-l) (3) 

for an electron wave packet ipn being initially a narrow Gaussian centered at an arbitrary lattice 
site Uq: 

ipn{0) = Aexp 

where A is the normalization constant. The solution of Eq. (3) can be expressed in terms of 
the eigenvalues X^, and eigenfunctions (p^n of the Hamiltonian (1) as follows 

N N 

Mt) = E E e"*^''V^^n<^^m^m(0) , (5) 
u=l m=l 

where the eigenfunctions (pi,n are chosen to be real. The quantities we use to characterize the 
dynamics of the electron wave packet are the mean position of the packet (centroid) : 

N 

= ^(n-no)lV'nWI' , (6a) 

n=l 

and its Fourier transform: 

1 poo 

f{u;) = — dte^^'xit). (6b) 

ZTT Jo 

3 Results and discussions 

In all simulations we used the helix of = 101 bases. The hopping integral was chosen to be 
J = 0.27 eV (see, e.g., Ref. [6]). The initial wave packet was placed in the middle of the hehx, 
no — 51. 

Figure 1 displays results of our calculations of the centroid x{t) for various magnitudes of the 
gating potential drop U±. At zero gating potential (upper panel), the motion of the centroid 
represents simple harmonic oscillations with the period tb — 27r/C/|| and the amplitude Lb ^ 
{4J/U\\)a, where 4J is the bandwidth, which is in full correspondence with the standard picture 
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Fig. 1. Bloch oscillations of the centroid x{t), Eq. (6a), calculated for a single-stranded helix of 101 

bases at various magnitudes of the gating field. The corresponding gating potentials are indicated in 
the plots. 

of Bloch oscillations in a hnear chain. The helical symmetry does not affect the oscillations at 
Ui =0. 



On turning on the gating field, U± ^ 0,, the motion of the centroid still manifests an oscillatory 
behavior, however, a more complicated one as compared to simple harmonic oscillations (Fig. 1 
middle and lower panels). The Fourier spectra / (a;) of the centroid plotted in Fig. 2 shed light on 
the situation. It demonstrates that a nonzero gating field gives rise to a splitting of the Bloch 
frequency into a multiplct with a frequency spacing dependent on the gating potential U±. 
The resulting signals presented in Fig. 1 (middle and lower panels) are formed because of the 
superposition of harmonic oscillations with several different and generally non-commensurate 
frequencies. 
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Fig. 2. Spectra of Bloch oscillations of the centriod depicted in Fig. 1. 
Finally, for the parameters we use, which are typical for the synthetic dry DNA, the period of 

oscillations = 27r/f/|| ~ 1 ps, i.e., it falls in the THz domain. 



4 Summary 

In summary, we have demonstrated that the intrinsic helix conformation of the DNA strands 
can have strong impact on its radiation properties. The electric field along the stack direction 
forces the injected charge to exhibit Bloch oscillations. In a tilted electric field, however, the 
harmonic Bloch oscillations become a superposition of oscillations with close and generally non- 
commensurate frequencies which can be tuned by an external electric field. The frequency of 
the oscillation falls in the THz domain. This finding is important for the self-assembled DNA 
arrays on gold with the DNA molecules being tilted with respect to the surface [13,14]. Such 
arrays may provide a nano-scaled source of coherent THz radiation. 
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